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Abstract 

We obtain the black body radiation formula of Planck by consid- 
ering independent contributions of multiphoton entities 
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1. Einstein introduced the light quantum hypothesis in 1905 ||T| by 
noticing that the entropy calculated from the Wien formula for the black 
body radiation 

I w (u, T) = av z exp(-bv/T), (1) 

which is the high-frequency, low-density limit of the complete Planck's for- 
mula 

,, ^ 8nhis 3 1 , . 

I (",T)= c3 e{hu/kT) _ v (2) 

this entropy had the same form as the entropy of a "gas" of independent 
particles, light quanta (later called photons). Viceversa, starting from a gas 
of independent light quanta one arrives only to Wien's formula, as shown 
since 1910 by many people ||. 

The question arises as how to get the complete radiation formula from 
the corpuscular concept of the photon. In this letter we present what we 
think is about the simplest derivation of Planck's radiation formula, starting 
from the idea of "multiphoton states" 0, i.e. molecules of light compound 
of n photons with zero binding energy. 

We know of course that the most satisfactory deduction of Planck's for- 
mula is obtained nowadays by incorporating the right B.E. statistics into the 
corpuscular picture, as first shown by S.N.Bose || and interpreted and ex- 
tended by Einstein in 1924/25 01]. Our calculation makes use of well-defined 
multiphoton states as the natural way of taking in account the identical, 
indistinguishable character of the photons. 

The idea of multiphoton "molecules" arose around 1910, after P. Debye 
presented |5j a simple, ondulatory method to obtain the radiation law; some 
historical comments are included below in this letter. 

2. We start from the photon entity as particles of energy e = hv (they 
are truly massless particles of helicity ±1 with energy and momentum e = 
hv = pc,p, but we want to think and write in the spirit of the 1910s). 

We make now this additional simple hypothesis : light of definite direc- 
tion and frequency v presents itself in units ("molecules") of 0, 1, 2, n, ...photons, 
with energy nhu, that is, with zero binding energy, and they contribute in- 
dependently to the energy density. 

The rationale for this assumption is, of course, in the spirit of the atom- 
istic or democritean point of view, that all the photons of same frequency 
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are created equal (identical), and that in an occamian, economical descrip- 
tion, we should abstain from considering situations in which the photons are 
individually distinct (indistinguishability); that is, there is a "state" with 
two photons different, for the counting, from the naive juxtaposition of two 
distinguishable photons; same thing for the three photon states, etc. So the 
multiphoton molecules should be thought of as entities only to the effect of 
counting, and not as claiming for the existence of some force which really 
binds the photons together. 

In this way one gets the correct radiation formula at once. We are unaware 
of any complete proof of this statement, so we now show the deduction, 
leaving some historical remarks for later. 

For the n-photon " molecules" of energy nhu, (n = 0, 1, 2, ...) the proba- 
bility is of course given by the Boltzmann factor 

p(n) = c exp(—nhis / kT) (3) 
^^£>(n) = 1 gives c = 1 — x, x = exp(—his/kT). (4) 

n=0 

The multiphoton states still have frequency is, so the common phase space 
factor is 2 (polarization) x47r(solid angle) xz/ 2 (radial 3D factor); hence the 
density of energy per inverse wavelength is 

I(is,T) = 87r(z/ 2 /c 3 )X~n/iz/(l - x)exp(-nhv/kT) 
8nis 3 hx 

c 3 1 — X 
Snis^h 1 
c 3 exp(hisfkT) - 1 K J 

i.e., the correct complete radiation formula. 

The two traditional limits of the radiation formula do have a particle in- 
terpretation, as already remarked by Wolfke [2]. E.g. at low temperature 
only the lightest "molecule" is excited, that is, the one-photon mode, of en- 
ergy his, that leads to the Wien's limit law (1), which was Einstein's starting 
point §. 

On the other hand, for high densities or temperatures all the n-photon 
states contribute together and sort of coalesce, giving a scale-invariant power 
law distribution I oc is 2 , i.e. the classical electromagnetic (Rayleigh-Jeans)limit 
formula. 
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3. If we now start from (5) and expand the denominator, we get 

oo 

I{v, T) = 8vr(z/7c 3 )(H exp{-nhu/kT) (6) 

n=l 

already obtained by several people in the 1910s; superficially, (6) looks like 
an addditive contribution from each n-photon state, and it was this analogy 
which prompted Ehrenfest |7|], Wolfke 0, de Broglie 0, |9| and others || 
to advance the "n-photon molecule" concept; notice however that the literal 
interpretation of (6) as a cooperative multiphoton formula falls short of the 
correct derivation in (5) on three accounts: 

1. There is no contribution from the zero-photon state 

2. The n— photon state contributes with energy hp, instead of the correct 
nhv 

3. Probabilities exp(— nhu/kT) do not add up to one. 

Still, it is remarkable that the expansions (5) and (6) are really identical! 

4. Why our simple derivation was not, to the best of our knowledge, 
obtained and exploited many decades ago, is for the historians to discuss; 
perhaps light quanta correlations due to indiscernibility was too novel at the 
time. Recall, however, that the Gibbs paradox was well known at the time, 
with the proposed solution by Gibbs himself [10]. Notice that this treatment 
of photons is totally correct even from the modern point of view; in particular 
the zero-point energy does not show up. Also we are really doing second quan- 
tization. Namely "first quantization" just amounts to counting frequency (or 
wavenumber) modes, because v = pc/h; the "second" quantization amounts 
to the existence of these photon molecules as separate entities. We finish by 
remarking that zero-energy bound states appear frequently in modern guise 
as Bogomol'nyi-Prasad-Sommerfield (BPS) states, and moreover uncharged 
PBS states, like the photon, are massless. 
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